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Abstract 

The Schrodinger functional (quantum/lattice field theory with Dirichlct boundary 
conditions) is a powerful tool in the non-perturbative improvement and for the study 
of other aspects of lattice QCD. Here we adapt it to improved gluon and quark 
actions, on isotropic as well as anisotropic lattices. Specifically, we describe the 
structure of the boundary layers, obtain the exact form of the classically improved 
gauge action, and outline the modifications necessary on the quantum level. The 
projector structure of Wilson- type quark actions determines which field components 
can be specified at the boundaries. We derive the form of 0(a) improved quark 
actions and describe how the coefficients can be tuned non-perturbatively. There is 
one coefficient to be tuned for an isotropic lattice, three in the anisotropic case. 
Our ultimate aim is the construction of actions that allow accurate simulations of 
all aspects of QCD on coarse lattices. 

1 Introduction 

It has become generally recognized in the last few years that accurate continuum 
extrapolations of lattice QCD will be possible in the foreseeable future only by using 
improved actions (see e.g. Q). 

In the Symanzik (on-shell) improvement program [||, ^, ||, [|, which we will follow, 
one systematically adds higher dimensional operators to an action, to cancel its lattice 
artifacts to some order in the lattice spacing a. Composite operators can be similarly 
improved. The coefficients of the improvement terms are trivial to calculate at tree-level; 
one-loop calculations are also possible, but already much more difficult, in general. Even 
if known, perturbative improvement coefficients per se are of limited value, since naive 
lattice perturbation theory does not work very well. It dramatically underestimates the 
improvement coefficients, for example. 

Tadpole improved perturbation theory J?]] works much better. For gluonic actions, 
where the leading errors of the Wilson plaquette action are 0(a 2 ), tree- level (or one- loop) 
plus tadpole improvement seems to cancel most of these errors, allowing quite accurate 



calculations on coarse lattices (o = 0.1 — 0.4 fm, say), as several studies || |T^, 11] 
have demonstrated. Similarly, it seems that the leading 0(a 2 ) errors of the Wilson or 
Sheikholeslami-Wohlert (SW) [0] quark actions, which break rotational symmetry, can be 
cancelled to a large extent with classical and tadpole improvement. This is demonstrated 
by the excellent dispersion relations [12, 13] of the D234 || actions, for example. 

On the other hand, for Wilson-type quark actions there are, at least at the quantum 
level, also 0(a) violations of chiral symmetry. Even without the explicit non-perturbative 
determination of the 0(a) improvement coefficient of the SW action (see below), there 
was clear evidence that this coefficient is not estimated sufficiently accurately by 
tadpole improvement on coarse lattices. In any case, it is obviously preferable to 
determine at least the leading improvement coefficients non-perturbatively. 

In seminal work, described in a series of papers (of which some are ] |l6| ||, |i~4{|), the 
ALPHA Collaboration has recently demonstrated how to do this. Using the Schrodinger 
functional [[l?], 18 and the demand that the PC AC relation hold for small quark masses, 
they determined the 0(a) improvement coefficient of the SW quark action (as well as 
other quantities) on quenched Wilson glue. They find that it is almost 20% larger than 
the (plaquette) tree-level tadpole estimate even for a relatively fine a = 0.1 fm lattice. 

Motivated by these observations, our aim in this paper is to adapt the Schrodinger 
functional to improved gluon and quark actions. 

The term "Schrodinger functional" is another name for the partition function of a 
euclidean quantum field theory with Dirichlet boundary conditions. In other words, 
one considers a theory with fixed boundary conditions on the fields at times xq = 
and xo = T, say. The Schrodinger functional is a powerful tool not just for the 
calculation of O(a) improvement coefficients of Wilson- type quark actions, but also for 
the determination of normalization and improvement coefficients of currents, as well 
as running couplings and quark masses. This is accomplished, in short, by using the 
finite volume and the boundary values of the gauge and quark fields as probes of the 
theory. Since Ward identities are local equations, they must hold independent of volume 
and boundary conditions. Probing the theory in this manner provides constraints that 
can be used to determine improvement coefficients. For the determination of running 
couplings one uses finite-size scaling techniques. One of the more technical reasons for the 
success of the Schrodinger functional is that the use of non-periodic boundary conditions 
in the time direction implies the absence of zero modes (at least for sufficiently small 
lattice spacings), allowing one to perform simulations directly at zero (or small) quark 
mass. This avoids a number of conceptual and practical complications. 

For the case of pure gauge theory in the continuum it is quite clear what is meant 
(at least in a formal sense) by Dirichlet boundary conditions, 

Z[C',C] = I IJA(x) J D[A]e- 3 °W, S g [A] = ^J d*x F^(x)F^(x) . (1.1) 

Here the path integral is over all gauge fields A fl (x), < xq < T, with boundary values 
(k = 1,2,3) 

, , , f C£(x) at x = 

A k ( x ) = { 7; ; (i.2) 

\ C k (x) at xq = T , 
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where, generally, C denotes the gauge field obtained after applying a gauge 
transformation A to C[] 

For lattice QCD it is not immediately obvious how Dirichlet boundary conditions 
should be defined. The Schrodinger functional for pure gauge theory, in the form of 
the Wilson plaquette action, was studied in [18]; Wilson quarks were treated in IS]. 
For the case of improved gluon and quark actions considered here, we will see that the 
"boundary" in the Schrodinger functional will consist of a double layer of time slices. 
This causes slight complications at intermediate stages, but ultimately the Schrodinger 
functional formalism for improved actions proves to be hardly more complicated than 
for Wilson gluons and quarks. 

There are several issues concerning the aims, scope and motivation of this paper that 
we should clarify at this point. First of all, the phrase "improved action" in the above is 
meant to refer to actions where at least the classical 0(a 2 ) errors have been eliminated 
by the addition of third and/or fourth order derivatives. Of course, in general it is 
inconsistent to consider improvement at order a 2 before having eliminated all errors at 
0(a). However, it is perfectly consistent to consider classical improvement to any order 
and then supplement it with quantum improvement to some lower order. It also makes 
sense to supplement classical with tadpole improvement .0 

For pure gauge theory perturbative improvement coefficients have been calculated to 
one loop for two actions [||, |2(J, so far. It is conceivable that at some point a practical 
method will be found to tune all gluonic a 2 errors to zero (at least for isotropic lattices) . 
For fermions, on the other hand, this is clearly impossible; there are just too many terms 
at order a 2 This leads to an obvious question: If we can only eliminate 0(a) but not 
all 0(a 2 ) errors, why should we even bother to eliminate what we hope are the leading 
a 2 errors via classical and tadpole improvement? We will argue that there are both 
theoretical and "empirical" reasons to suggest that this is a worthwhile enterprise. Let 
us discuss these reasons in turn. [The reader not interested in motivational issues, but 
only in the Schrodinger functional per se, can at this point skip ahead to the outline at 
the end of this section.] 

Remember the special nature of the 0(a) errors of Wilson- type quarks. They arise 
because we added a second-order derivative term to solve the doubler problem. This 
term breaks chiral symmetry, and, at least at the quantum level, 0(a) errors become 
unavoidable.^] On the other hand, these terms do not break rotational symmetry at 0(a) 
or 0(a 2 ). This is in contrast to the leading a 2 errors of both gluon and quark actions, 
which break rotational but not chiral symmetry. One might therefore argue that it is 
somewhat besides the point to label these errors as 0(a) and 0(a 2 ) errors. Viewing them 

lr The integral over gauge transformations A(x) at xo = (or xq = T) is necessary to ensure the gauge 
invariance of the Schrodinger functional. 

2 Note, though, that for dynamical QCD, where the quark dynamics feeds back into the glue, it might 
be necessary to use an iterative procedure to self-consistently tune the tadpole factors and the O(a) 
coefficient (s) of the quark action. It seems likely that this iteration will have converged to sufficient 
accuracy after just a few steps. 

3 If other terms are suitably added by a field transformation the chiral symmetry breaking can be 
avoided at the classical level to basically arbitrarily high order; cf. [ft and appendix C. 
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instead as the leading violations of chiral, respectively, rotational symmetry, it makes 
perfect sense to try to eliminate both of them as best as one can. 

Now to the empirical reasons. Here the basic point is simply that the elimination 
of the leading 0(a 2 ) errors through classical and tadpole improvement works quite well. 
We already referred to the surprisingly accurate results obtained with improved gluon 
actions on coarse lattices. Concerning the more difficult case of 0(a 2 ) improved quark 
actions, I argued in [15], after surveying the available data, that all quantities with a 
weak dependence on the 0(a) improvement term (dispersion relations, mass ratios) are 
also accurate to a couple percent on coarse lattices when calculated with such an action. 
It is only quantities with a relatively strong dependence on this term (like the p mass 
in units of a gluonic scale such as the string tension, or the hyper-fine splitting of heavy 
quark systems) that still have significantly larger errors on the coarsest lattices. It is 
therefore not unreasonable to believe that with no 0(a) and only small 0(a 2 ) errors 
accurate simulations of QCD will be possible on coarse lattices. 

Concerning the gluon actions there is another more technical, but potentially 
important, issue we should mention. Namely, in addition to the 0(a 2 ) errors of the 
Wilson plaquette action, there also seem to be large non-power-like lattice artifacts 
that set in on lattices with spacings slightly above 0.1 fm. What we are referring to 
is the non-trivial phase structure of the plaquette action when embedded in a larger 
space of couplings (fundamental and adjoint representations of the plaquette). This 
was the subject of many studies in the early days of lattice gauge theory; cf. [21] for 
recent results and references. The (in)famous dip of the discrete /^-function, a similar 
dip in the ++ glueball mass, and non-smooth behavior in various thermodynamic 
quantities are all believed to be consequences of a second-order phase transition close 
to the fundamental coupling axis of the Wilson gauge action. Although further studies 
are needed, there is evidence that these lattice artifacts are significantly smaller with 
improved gluon actions (see e.g. [|ll|, [2^]). There are also more direct indications that the 
second-order phase transition point moves away from the fundamental coupling axis with 
improvement [23]. Finally, it seems that the variance of the static potential measured 
with improved glue §| |24| is much smaller than for Wilson glue at the same lattice 
spacing. This might also be related to the phase structure (on the other hand, there 
might be more prosaic reasons for this; this remains to be investigated). 

The potential significance of these issues is as follows. In their study of the 0(a) 
improvement of the SW action on quenched Wilson glue Liischer et al found [14] that 
the improvement coefficient could not be determined for lattices coarser than 0.1 fm, due 
to large fluctuations ("exceptional configurations") of the glue, leading to accidental zero 
modes at small quark mass. In the quenched approximation this problem is bound to 
occur at some sufficiently coarse lattice spacing, but it is not clear a priori exactly where 
it should occur. It could depend on the specific lattice action used. The observations 
described in the previous paragraph, lead one to hope that with improved glue^| the 

4 And perhaps order a 2 improved quark actions, which in our experience have less problems with 
"exceptional configurations". (Which, by the way, means that the latter term is strictly speaking a 
misnomer; the occurrence of near zero modes depends also on the specific quark action used for a given 
configuration.) 
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breakdown will occur on coarser lattices. 

While finishing the write-up of this paper, hard facts began to emerge that support 
the above arguments: As will be described in [25], the breakdown of the quenched 
approximation does indeed occur on coarser lattices when improved glue is used. 

The hope that this would happen was the main motivation underlying this paper. 
To obtain accurate continuum results from lattice QCD requires simulations at a series 
of lattice spacings that provide a significant "lever arm" for extrapolations. Due to the 
dramatic increase in the cost of a QCD simulation as the lattice spacing is decreased, 
such a series can be obtained orders of magnitude cheaper if one makes use of coarse 
lattices. 

Concerning the scope of this paper we should mention that we will also consider 
anisotropic lattices || |26|, 27]. Such lattices, with a smaller temporal than spatial lattice 
spacing, offer the best hope to tackle certain hard problems, like the simulation of heavy 
quarks in a relativistic formalism |b|, || and detailed studies of glueballs ]p^| . Due to 
the breaking of manifest space-time exchange symmetry on anisotropic lattices, there 
are more terms that have to be tuned to eliminate all 0(a) errors of a Wilson- type quark 
action (as we will discuss). Tuning these errors to zero is therefore more complicated, 
but should be within reach of the Schrddinger functional technology (cf. sect. ^). Note 
that anisotropic lattices fit quite naturally into the Schrddinger functional framework, 
since the boundary conditions already treat space and time differently. 

The outline of this paper is as follows. In sect. 2 we derive the Schrddinger functional 
for improved gluon actions. We discuss the structure of the boundary layers and the 
precise form of the (classical) action at the boundaries, where the fields are given fixed 
values. We then specialize to spatially constant boundary values and discuss the classical 
background field they induce in the bulk, including its uniqueness. The subject of 
sect. 3 is Dirichlet boundary conditions for (improved) quark actions. All Wilson-type 
quark actions of practical interest, on isotropic as well as anisotropic lattices, have a 
certain projector structure, which dictates which field components can be fixed at the 
boundaries. Once these are specified, one can define boundary quark fields as functional 
derivatives with respect to the boundary values. These boundary fields are extremely 
useful in applications of the Schrddinger functional. One application is described in 
sect. 4, namely, the determination of the 0(a) improvement coefficients of Wilson-type 
quark actions. Concluding remarks and a brief discussion of future work can be found 
in sect. 5. 

Appendix A summarizes some of our notation and conventions. In appendix B 
we discuss the classical boundary errors of the pure gauge Schrddinger functional. In 
appendix C, finally, we derive the form of 0(a) on-shell improved quark actions on 
anisotropic lattices. We also recall |(| the classical improvement of quark actions. 

This paper is largely self-contained, but the reader might want to consult fl8[ for a 



detailed introduction to the Schrddinger functional, ref. [16] for a summary of various 



applications, and p, n] for discussions of Symanzik improvement for fermions. 
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2 Gauge Fields 



2.1 Gauge Actions in the Bulk 

Consider a four-dimensional hypercubic lattice of extent and lattice spacing in 
direction /j, = 0, 1,2,3. Except for purely classical considerations we assume the lattice 
to be spatially isotropic with a = a^, k = 1,2,3. Other notations and conventions are 
summarized in appendix A. 

To avoid having to worry about boundary effects at this point, let us impose periodic 
boundary conditions on the gauge fields for the moment. The lattice gauge actions of 
interest can be written as 

S M = J;E E c r ReTr[l - [/(C)], (2.1) 

where U(C) denotes the product of link fields along the loop C, r labels a finite set 
of types of loops, and S T is the set of all positively oriented (say) loops of a given 
type that can be drawn on the lattice. On an isotropic lattice the types of loops might 
include plaquettes, rectangles, etc; on an anisotropic lattice one would have to distinguish 
spatial plaquettes, temporal plaquettes, spatial rectangles, short-temporal rectangles, 
long-temporal rectangles, etc. 

The real numbers c T are determined by some improvement condition. We do not 
really need to know the exact numbers to set up the Schrodinger functional, but for 
concreteness let us be more explicit about a simple class of actions that is sufficient to 
illustrate what happens for more general improved actions. Namely, for "plaquette + 



rectangle" actions [28] we can write 



(0) 1 - PjjAx) _ (1) 1 - R^{x) _ (1) l-R Ufl (x) 



W = 3 E ^ ReTV 

(2.2) 

where AV = JT„ a M is the volume of a single lattice hypercube, and and R^ u are 
equal to U(C) with C a 1 X 1 (plaquette), respectively, 2x1 (rectangle) Wilson loop in 
the [iv plane. The coefficients cj2) = 1 + 4(c~^ +c„ ) and for good measure we have also 
indicated tadpole improvement factors that one might want to include^] As special 
cases this class of actions contains the Wilson and tree- level 0(a 2 ) improved [28, [J 
actions, 

,„ f Wilson action 
c { u = { i (2-3) 
I 12 Improved action . 

Another interesting case, on anisotropic lattices, is defined by Cp = jk(l — 5^o) (cf. p7|). 
Only the spatial directions are improved now; there are no tall-temporal rectangles in the 
action, leading to Oq errors. Since this action contains only loops of maximal temporal 
extent 1, there are no unphysical branches in its dispersion relation. 



On anisotropic lattices additional renormalization factors might have to be included, depending on 
how one defines the tadpole factors. 
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2.2 Transfer Matrix 



The Schrodinger functional can either be defined as the partition function with fixed 
boundary conditions, or as a suitable power of the transfer matrix evaluated between 
initial and final states. On the lattice the meaning of "fixed boundary conditions" is not 
a priori clear. Understanding this issue via the transfer matrix is at least pedagogically 
useful, so we will proceed by reviewing the transfer matrix formalism. For improved 
gauge actions this formalism was described in p9| ], whose approach we will follow. 

We will use the plaquette + rectangle actions as examples, but as we proceed it will 
become clear how to construct the transfer matrix for more general actions. We will also 
see that the general form of the transfer matrix is determined solely by the maximal 
temporal extent of the loops appearing in the action. To avoid tedious and irrelevant 
case distinctions, we will use the term "Wilson action" as a shorthand for any gauge 
action with loops of maximal temporal extent 1. Note that this also includes plaquette 
+ rectangle actions without tall-temporal rectangles, which, as mentioned earlier, are 
useful for anisotropic lattices. We will use the phrase "improved action" to denote actions 
with loops of maximal temporal extent 2. 

For pedagogic reasons we first consider the Wilson case. The classical lattice field 
equations for such actions (in temporal gauge Uq{x) = 1, say) are second order difference 
equations in xq. Specifying the link field on one initial and one final time slice leads to a 
unique solution to the field equations. The transfer matrix of the quantum theory should 
therefore give the amplitude to go from some gauge field on one time slice to another 
field on some other time slice. 

To define the transfer matrix we first have to specify the Hilbert space Ti on which 
it is supposed to act. We take Ti to consist of wave functions <f[U] that are functions of 
the gauge field U on a fixed time slice, with the natural inner product 

(H^) = J D[U] <p 2 [U\* tp x [U] , D[U] = [] dU k (x) , (2.4) 

in terms of the Haar measure dU on SU(iV). More precisely, Ti consists of all gauge 
invariant wave functions of finite norm. Note that any wave function can be made gauge 
invariant by averaging over gauge transformations U — ► U A . This defines a projection 
operator onto the physical subspace, 

P<P[U\ = / II dA (xM£/ A ] • (2-5) 

J X 

We can write (p[U] = (U\ip) where the states \U) are eigenstates of the link field operators 
Uk(x) in the Schrodinger representation, with eigenvalues Uf~(x). For gauge invariance 
the projector IP should be applied to the |J7)'s, which will be understood in the following. 
The identity can be decomposed as 1 = J D[U] \U)(U\. 
The transfer matrix T acts on wave functions as 

(J<p)\U'] = J D[U] (U'\T\U)<p[U\ , (2.6) 
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Figure 1: Temporal weight factors of various loops for Wilson (left) and improved 
(right) action transfer matrices. 



where its matrix element (or kernel) to go from gauge field U at time t to U' at time 
t + oo is given by 

(U'\T\U) = f l[dV (t,x) e - AS V> v °M , Vo (t }X ) G sU(iV) . (2.7) 

Here AS[U', Vq, U] is given by exactly the same expression as the action itself, restricted 
to a two time slice lattice with gauge fields U at t, U' at t + «0i an d Vq playing the role 
of the temporal link fields^ — except that the contribution from loops completely within 
the spatial planes at t and t + a® are weighted with a factor of ^. This is illustrated in 
figure |[ It is clear why this factor of | must be present: When we take powers of the 
transfer matrix (to get the partition function or Schrodinger functional, see below) all 
gauge fields on the internal time slices contribute twice. Similarly, a correlation function 
involving gauge fields is equal to some power of the transfer matrix with operators U 
inserted at appropriate places in the product. 

We now turn to the 0(a 2 ) improved case. With loops in the action extending up to 
two steps in the temporal direction, the action involves second order, the field equations 
fourth order temporal differences (in temporal gauge). There will be a unique solution 
to the field equations if we specify the gauge field on two initial and two final time 
slices. Correspondingly, we expect the Hilbert space to consist of wave functions on two 
consecutive time slices, and the transfer matrix to give the amplitude for transitions 
between such double layers. 

Let 

<p[U] = ^[U^\U (0 \U^} (2.8) 

be a wave function depending on the link field on time slice t, on time slice 
t + do, as well as the temporal links connecting the two time slices. Let the Hilbert 
space consist of gauge invariant wave functions of finite norm, with respect to the inner 



6 Note that the integral over Vo will project ip in (2.6) onto its gauge invariant part, in case it was not 
gauge invariant to start with. 
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product 



M<Pi) = / D[U] <p 2 [U\*<Pi[U] , ^] = n[ dC/(0) ( x )Il dul + \x)dul-\x) . 

(2.9) 

The transfer matrix is then defined by 

(U'\T\U) = e-^ u 'M J] • ( 2 -10) 

x,fc 

The 5-function (with respect to the Haar measure) means that the transfer matrix is 
effectively a "2 — ► 1" and not a "2 — > 2" time slice operator. AS[U', U] is the action of a 
three time slice lattice obtained by fusing the top layer £A + ) f U with the bottom layer 
[/'(-) of £7', where we weigh a Wilson loop of temporal extent At (in lattice units) by a 
factor 

. ' 6 at xo = t 

w(At = 2) = l, w(At = l) = -, w(At = 0) = < 1-26 at x = t + a (2.11) 

I 6 at xo = t + 2ao 

This is presented graphically in figure [I]. Note that {U'\T\U) is gauge invariant, since 
AS is. 

Time reflection symmetry dictates the above form of the "temporal weight factors" 
w(At). When taking powers of the transfer matrix the w(At) combine to give 1 for all 
loops in the bulk, ensuring that we obtain the correct (bulk) partition and correlation 
functions. This is true for any value of the (real) parameter b. The effect of b is only 
noticeable at the boundaries. Its value is fixed, classically, by demanding the absence of 
order a 2 boundary errors.[] This will be discussed in the next subsection and appendix B. 
Let us for the moment keep b a free parameter. 

The transfer matrix of improved actions is not hermitean J29| ], since energy 
eigenvalues are complex at the order of the cutoff (this can also be seen by studying the 
dispersion relation). In most situations the resulting effects, chiefly oscillatory behavior 
of correlation functions at small times, do not create significant problems in practice]^] 

2.3 Gauge Actions with Dirichlet Boundary Conditions 

We can now write down gauge actions with Dirichlet boundary conditions. We choose the 
spatial boundary conditions on the gauge fields to be periodic. To obtain the Schrodinger 
functional on a lattice of temporal extent T/an we just have to evaluate T T//a ° between 
fixed initial and final states. In the Wilson case the gauge fields have to be specified at 
xq = and xq = T, in the improved case on the xq = — ao, and xq = T',T' + ao double 

7 In periodic boundary conditions were considered, where the issue of boundary errors does not 
arise, and the "symmetric" choice b = 1/3 was used. As we will see, this is not the value minimizing 
boundary errors with Dirichlet boundary conditions. 

8 The one obvious exception being glueballs, whose bad signal to noise properties force one to work 
with the first few time slices. As mentioned in the introduction, anisotropic lattices can help here. 
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Figure 2: Temporal weight factors of various loops in the Wilson (left) and improved 
(right) actions with Dirichlet boundary conditions on lattices with L = 5a and T = 5ao- 
Time, measured in lattice units, runs upwards, as marked on the left (T 1 = T — 1 for the 
improved case). The thick lines indicate the boundaries of the bulk region; every loop 
with at least one link strictly inside the bulk has a temporal weight factor of 1. 



layers. To allow for a unified notation for the Wilson and improved cases in various 
formulas, we have introduced 

f T Wilson action 

T = { (2.12) 
[ T — ao Improved action . 

Denoting the boundary gauge fields by W and W', we define the Schrodinger functional 
for pure gauge theory as 

Z[W',W] = {W'\T Tla °\W) . (2.13) 

Remember that in the improved case W is a shorthand for the fields W^ + \ W"H 
on a double layer (and ditto for W'). More explicitly, we can write 



Z[W',W] = J D[U) e S,[U\W',W\ ) D[u] - -Q 



T'-a 

du (o,x) n Hdu^x) 

xo=ao 



S g [U;W',W] = E w(C)crBBTt[l - U(C)], (2.14) 

9 r C£S T 

where the temporal weight factor w{C) for a loop in an action with Dirichlet boundary 
conditions is best defined pictorially, as in figure |2[ 

As remarked in the caption of this figure, all loops containing a dynamical link have 
a temporal weight factor of 1. This fact simplifies the coding of the updating algorithm 
for the Schrodinger functional. 
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We now discuss the boundary errors of the Schrodinger functional. In particular we 
have to fix the so far free parameter b in the improved case. In view of the remark in the 
previous paragraph, the reader might wonder at this point, how errors associated with 
the boundary fields, which are, after all, fixed in a simulation, can have any effect on an 
observable defined in terms of dynamical variables in the bulk. 

First of all, even on the classical level it is legitimate to consider observables defined 
by taking functional derivatives with respect to the boundary values, and only then 
fixing them at some suitable value. This introduces cross terms between the bulk and 
the boundary fields. The "SF coupling" defined in p0[] , for example, does involve such 
derivatives. Secondly, on the quantum level the boundary effects can "spill into the bulk" 
and affect the dynamics there, as we will see shortly. 

Let us first consider the classical boundary errors of the Schrodinger functional. 
We want to know the difference between the lattice and the continuum actions with 
corresponding boundary conditions. This difference has two sources. One is the 
difference between the continuum and lattice action densities (the integrand/summand), 
the other the difference between an integral and a discrete sum approximating it. 

This is basically a classical calculus problem; the elementary, if slightly wordy, 
tracking of the differences for the Wilson and improved actions is discussed in appendix B. 
We here only state the conclusion: The Wilson action has order a 2 boundary errors, and 
the improved action only a 4 errors if we choose b = 1/24. 

Actually, at tree level the precise value of b does not matter for the special (spatially 
constant) boundary values of the gauge field we will consider later. For many applications 
this is even true on the quantum level at 0(a); a point to which we will return in a 
moment. 

In the quantum case it is expected [17, 18] that in addition to the usual counter terms 
in the bulk, some small number of additional counter terms might be necessary on the 
boundaries. In pure gauge theory the only possible terms at 0(a) are spatial sums over 
the boundary time slices of lattice versions of TrF 2 fc and Tri^. Note that there is no 
difference in this respect between isotropic and anisotropic lattices (in contrast to the 
bulk, where fewer counter terms are necessary in the isotropic case). 

So, no terms not already present in the classical action are required on the quantum 
level. In addition to the usual renormalization of the gauge coupling in the bulk (and 
other bulk renormalizations for anisotropic lattices; see 27], and for first one-loop 
results [31]), there will be corrections to the weight factors of some of the loops at 
the boundary. It is interesting to note that to obtain the counter term for Tr F^ k in 
the Wilson case by necessity affects the temporal weight factor of plaquettes involving 
dynamical links (namely the temporal plaquettes touching the boundaries at xo = and 
xq = T) , whereas for the improved case one can choose to renormalize only the weight 
factors of temporal plaquettes within the boundary doubler layers, which do not involve 
dynamical links. To obtain the counter term for Tri 7 ^ one can, in both cases, suitably 
adjust the weight factor of spatial plaquettes within the boundary.^] 

This means that for the improved case some observables will not be sensitive to 



3 For the currently known one-loop corrections to the boundary weights cf. j^] and references therein. 
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0(a) quantum errors incurred by the wrong choice of the temporal and spatial plaquette 
coefficients in the boundary; only observables involving derivatives with respect to the 
boundary values that lead to a coupling of the dynamical links with the boundary 
plaquettes will "see" such errors. 

In any case, though, for most applications of the Schrodinger functional it is not 
important to know the boundary counter terms. The basic reason has already been 
mentioned in the introduction: The boundary fields are used to probe the local dynamics 
of the theory, which, to leading order, should be independent of the precise details of the 
boundary terms, including their renormalization (see || for a formal proof). 



2.4 The Classical Background Field 

The boundary values of the gauge field induce some classical background field in the bulk. 
It is important, e.g. for perturbative calculations, that the background field be unique 
(up to gauge equivalence). This will not be the case for an arbitrary choice of boundary 
values. Spatially constant boundary values have turned out to be both convenient and 
sufficient for all applications of the Schrodinger functional (so far, at least). Physically, 
these boundary conditions correspond to a constant color-electric background field. We 
will first consider generic boundary conditions of this type and later discuss further 
restrictions that have to be imposed when we turn to the question of uniqueness. 

In the continuum the background gauge field in such a situation, denoted by A^(x), 
is given by 

A (x) = , A k (x) = i [x C' k + (T- x )C k ] . (2.15) 
Here C k and C' k are the (hermitean) boundary values of the gauge field, 

t / \ \ C k at x = 

M*) = \ r , . „ (2-16 
[ C k at x = T , 

in terms of which the color-electric field strength reads FQ k = (CI — C k )/T. 

If one assume that C k and C' k commute (e.g. if they are diagonal, see below), then 
an obvious guess for the background field on the lattice is U^(x) = exp(— ia^A^x)), or 
explicitly, 

U (x) = 1 , U k (x) = exp[-i^[x C' k + (T — x )C k ]) . (2.17) 

For the Wilson case this corresponds to the boundary values 

f W fc (x) = e - iakCk at x = 

Udx) = { V ' . r , (2.18) 

V ' \ W^x) = e- ia * c * at x = T . 17 



For improved actions one can read off the boundary fields wjf 1 ^ and from (|2,1 

at t = -o ,0 and t = T',T' + o , respectively. W^{x) = 1 = W'^(x), of course.f" 

10 Note that we are reverse-engineering the boundary conditions: We start with a gauge field that, as 
we will see, satisfies the lattice field equations, and read off the boundary conditions from the boundary 
values of the field. This is perfectly legitimate. 
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It remains to be proved that ( 2.17 ) actually is the correct background field on the 
lattice. To this end, note that for this field the non-vanishing contributions to the action 
come from the temporal plaquette and rectangle terms, 

P ok (x) = P k0 (x)* = e U-i^[C' k -C k ] 



T 

R k(x) = R k0 (x)* = exp(-2^K-C fc ]) . (2.19) 

The contribution of a loop to an arbitrary gauge action is simply related to the color- 
electric flux through the loop. Since the background electric field is constant in space 
and time, the same is true for these fluxes, which implies that ( p. 17 ) obeys the lattice 



field equations for any gauge action. The uniqueness of the background field ( 2.17j ) will 
be discussed in the next subsection. 

The fact that ( |2.17| ) satisfies the field equations for an arbitrary gauge action, testifies 
to the "quality" of the abelian boundary fields, as previously remarked in [O]. This can 
also be seen by calculating the value of the plaquette + rectangle action ( |2.2| ) for the 
background field with Schrodinger functional boundary conditions, 



S 9 [U] = ^^ETrf^-sina^fl + ^+ci^sin 2 ^) 

= ^^E^-C^ + OCa 4 ), ak = ^(C' k -C k ). (2.20) 

" k 

The same result would have been obtained with periodic boundary conditions in the 
time direction. We see that for constant abelian boundary fields even the Wilson action 
has only 0(a 4 ) errors relative to the continuum limit. 



2.5 Uniqueness of the Classical Background Field 



It is important to know if, for given boundary values of the gauge field, the background 
field is the unique (up to gauge equivalence) absolute minimum of the action. 



In ref. [18] a theorem establishing uniqueness was provided for the Wilson plaquette 
action, in any number of dimensions, for a certain class of diagonal and spatially constant 
boundary values. To specify this class consider diagonal boundary values 



Cfc 



1 

Lk 



I 4>kl 

o 



V o 





4>k2 



■ \ 

o 

• 4>kN I 



(2.21) 



and similarly for C' k in terms of a set of 4>' ka . Gauge symmetry implies that the cj) ka 
and (fi' ka have angular character; we can always choose them to lie in (— tt, it]. A vector 
(4>i, . . . , 4>n) is said to lie in the fundamental domain if 



E 



o, 



>i < <h < 



<4>n , 4> 



N 



< 2ir. 



(2.22) 
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It was proved in [18|] that a sufficient condition for the uniqueness of the background 
field fl2.17|) is that the vectors 4>k, 4>' k (defined in the obvious way) lie in the fundamental 
domain.^] The theorem also holds if some or all of the <f>k or (j)' k are zero (vectors). 

The proof given in [^] for the Wilson action is quite involved. As a corollary it 
implies that the analogous theorem holds in the continuum. One therefore expects such 
a theorem to hold for improved actions, since they are supposedly more continuumdike 
than the Wilson action. However, there is at least one step in the proof, related to the 
use of the field equations, that prevents a straightforward extension to improved actions. 

Even though uniqueness presumably holds for all "reasonable" improved actions, we 
have no proof at the moment. We have therefore investigated the question of uniqueness 
numerically. For various actions, volumes and boundary conditions we started with 
some random initial configurations, and used simulated annealing to find a minimum 
of the action. Besides the plaquette + rectangle action, we considered the action 
with an additional "parallelogram" term, i.e. a set if loops sufficient for full order a 2 
improvement ||. We stepped through the region of improvement coefficients of interest 
in practice — and beyond. Several boundary values in (and on the boundary of) 
the fundamental domain were considered. In all cases where the coefficients of the 
improvement terms are such that positivity holds (cf. Q; any viable action must satisfy 
positivity), the simulated annealing algorithm eventually converged to the classical 
background field Q2.17| ). 

We should add that we found the above to be true even on the very small lattices, 
like L = 2a, T = 3a,0 where the proof of [lq] does not apply. The only cases we could 
find where the global minimum is definitely not given by ( 2.17 ) involve actions with 
large negative coefficients for the improvement terms, violating positivity. Then indeed 
rather interesting things can happen. The simplest, and surprisingly rich, example is 
provided by the plaquette + rectangle action on a two-dimensional lattice with L = 2a, 
T = 3a and vanishing boundary fields, = = 0'. For a certain region of large negative 
rectangle coefficient, for instance, there is a discrete set of degenerate vacua that break 
spatial translation symmetry. The vacua have an anti- ferromagnetic character, in that 
translation invariance holds only under shifts by two lattice units. 

We conclude that for boundary values in the fundamental domain all improved actions 
of interest in practice seem to have the expected unique global minimum for arbitrary 
volumes. 



11 Another technical requirement in the proof is that the space-time volume is not too small. The lower 
bound given in reads TL > 20ao<2 for SU(3) gauge fields on a spatially symmetric lattice. 

12 The case T = 2a is trivial for improved actions, since there are no dynamical spatial links then. But 
we checked T — 2a for the Wilson case. 
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3 Fermions 



3.1 Fermion Actions in the Bulk 

Before writing down the most general quark action 

S g [ijj,ip,U] = a a 3 ^2 tp(x)Qip(a 



(3.1) 



we will be interested in, let us consider two important special cases on isotropic lattices. 
These are the "isotropic SW" M and the "isotropic D234" || actions, whose quark 
matrices Q can be written as 



Qwi 

Qw 2 



m 
m 



in terms of the "Wilson operator" 



(3.2) 



(3.3) 



We will therefore also refer to these actions as the Wi and W2 actions, respectively, as 
anticipated in eqs. (3.2). 

In the above a ■ F is the so-called clover term (cf. appendix A), and we used the 
standard first and second order covariant derivative operators, defined, on a general 
lattice, as0 



2a,, 



U,j,{x)if>(x + (J,) - U-,j,(x)ip(x - fi) 
U^{x)^{x + n) + U^(x)ip(x - y) - 2il>(x) 



(3.4) 
(3.5) 



(cf. appendix A for unexplained notation). The main point of the Wilson operator is 
that it can be expressed in terms of two-dimensional (spinor-) projection operators as 



P± = 1(1±7m) 



where we introduced forward and backward covariant derivatives, 



Vt^(x) 



l 



U±^(x)ip(x ±n) - ip(x] 



(3.6) 



(3.7) 



This property allows for an efficient coding of the Wi and W2 actions. 

On the classical level, where the clover coefficient u = 1, the Wi action has a 2 , the 
W2 action a 3 errors, and they are the "cheapest" isotropic actions with these errors. 



13 We ignore tadpole factors for the link fields that one might want to use in the context of an actual 
simulation. 
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They are natural candidates for 0(a) quantum improvement, which involves tuning the 
clover coefficient uo non-perturbatively. 

As we will see, to be able to define simple Dirichlet boundary conditions all we need 
is that the temporal parts of the quark actions have the above projection properties. On 
anisotropic lattices that are spatially isotropic, the actions of interest can therefore be 
written as 

Q = 7o W + boaoWo 2 + cB + M. (3.8) 

We have separated the terms in the quark matrix according to their 7-matrix structure. 
For the spatial part of ( j3ij| ) we need to make no special assumptions, so D = J2k JkDk 
is written in terms of generic derivative operators D k = + 0(a^)J^ and M denotes 
all other terms, 

M = m - T -^- ( iv ^2 <?Q k F ok + lo °kiFki + J! A fc + . . . J . (3.9) 

^ k k<l k ' 

The natural generalization of the Wi and W2 actions to anisotropic lattices, which we 
will also refer to by these names, involves 

Wi action: b§ = , r = 1 

W 2 action: b = - , r = - , (3.10) 
o 3 

as on isotropic lattices. 

As shown in appendix C, quantum 0(a) on-shell improvement of such actions requires 
the tuning of three coefficients: The temporal and spatial clover coefficients, loq and u, 
as well as c, the "bare velocity of light" , that has to be adjusted to obtain a renormalized 
velocity of light of l.P] 

We should add a remark on notation. For the isotropic Wi, or SW, action csw = ru 
is often referred to as the clover coefficient. It should not cause confusion when, more 
generally, we use that name for oj. 



3.2 Field Equations and Dirichlet Boundary Conditions 

Dirichlet boundary conditions for Wilson quarks were derived in |l^] using the transfer 
matrix formalism |32| (cf. also [|33|| ). This is somewhat involved, especially for improved 
quark actions. Since ultimately we do not need the transfer matrix itself, we will use a 
shortcut. Namely, a little bit of thought shows that to decide which field components can 
be consistently fixed at the boundary it is sufficient to consider the projector structure 
of the quark action (or field equations). 

14 Usually we use D M to denote the covariant continuum derivative, but its use in this section as a 
generic lattice derivative should not cause any confusion. 

15 In the case of dynamical simulations on anisotropic lattices one might have to iteratively tune 
parameters in the glue and the quark actions until they are self-consistent, i.e. correspond to the same 
renormalized anisotropy. 
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Figure 3: Dependence diagrams of P + ip (left side of each figure) and P~ip (right side 
of each figure) for Wi (left) and W2 (right) equations of motion. The field components 
pointed to by the arrows are determined by the field components in the time slice(s) at 
which the arrows originate. The corresponding diagrams for the ipP^ components are 
obtained by reflection with respect to a vertical line through the middle of each diagram 
(or by exchanging + and — at the top). 



The lattice equation of motion of ip in the bulk is derived by varying the action 
defined by ( [SO] ) with respect to ?/>. Projected onto P^ = P^ 1 components the field 
equations read 

P±il>(x) = -M ± P ± ^{x) - cJ2lkD k pTiP( x ) , < x < T' , 

k 

(3.11) 

with M ± = P ± MP ± . Analo gous equations are obtained for the ip components. 
Diagrammatically the structure of these equations is illustrated in figure ||. When one 
now asks how to fix some field components at the boundary in a way that is consistent 
- leading to a system of equations that is neither under- nor overdetermined — one 
finds the situation depicted in figure ||.[^] 

This implies that for a Wi action consistent Dirichlet boundary conditions correspond 

16 In figures ^ and [| we assume a fixed gauge field background. This is compatible with the standard 
method of performing the path integral in lattice gauge theory: first integrate out the fermions 
analytically for fixed gauge field, then perform the path integral over the latter by Monte Carlo. We 
mention this because otherwise the (standard form of the) clover term hidden in M would lead to 
complications; it connects to time slices beyond the ones shown in figure []. Also, this is another reason 
to avoid using the transfer matrix for a combined fermion and gauge system. Due to the clover term it 
would in general be more complicated than simply the (symmetrized) product of the gauge and fermion 
transfer matrices. This lead the authors of jw) to propose a modified clover term that extends only over 
two time slices instead of three. 



±V^ + 6 a (V^) 2 
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Figure 4: Schematic representation of the fermion field equations with Dirichlet boundary 
conditions for a Wi (left) and W2 (right) action. For clarity the field equations are 
represented only by the "diagonal" arrow of figure ||; what matters is that one arrow 
corresponds to one equation for a field component in terms of field components on other 
time slices (represented more correctly by two, respectively, three arrows in figure ||). 
This implies that the a priori unknown field values in the bulk (open circles) are uniquely 
determined by the boundary values (solid circles) and the field equations; there are 
exactly as many unknowns as equations. The corresponding diagrams for the 
components are obtained as in figure ||[ 
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to specifying 



P + i/j{x) = 


P(x) 


at 


XQ 


= 


${x)p- = 


P(x) 


at 


Xq 


= 


P~ip{x) = 


P'(x) 


at 


XQ 


= r 


^{x)P + = 


P'(x) 


at 


Xq 


= r 



(3.12) 



as derived in [19] using transfer matrix methods. Note that it makes sense to couple Wi 
quarks to either Wilson or improved glue. 

W2 quarks, on the other hand, can be consistently coupled only to improved glue 
with Dirichlet boundary conditions; this will always be understood in the following. In 
this case we have to specify not only the field components indicated in ( |3.12 ) on the 



inner boundary layers, but the same components also on the outer boundary layers at 
xq = —cio and T' + ao (cf. figure §). To avoid cumbersome notation we will not introduce 
additional symbols for the boundary values on the outer boundary layers. It will always 
be clear from the context if p, p etc refer to the inner or outer boundary, or both. 

All components at the boundaries not of the general form shown in eq. dUp have 



no influence on the bulk dynamics and could in principle be set to arbitrary values. 
However, for notational convenience we will set them to 0, and formally define ip(x) and 
tp{x) for all —00 < xo < 00 by "padding with zeros". That is, we set ijj(x) and ip(x) to 
zero for xq < and xo > T in case of a Wi action, and for xq < —oq and xq > T' + ao 
in case of a W2 action. This has the advantage that we can now define the lattice quark 
actions with Dirichlet boundary conditions to be given by exactly the same sum as for 
T = 00, with the boundary values as discussed above. 

Actually, we have cheated here. All we have shown above is which field components 
are to be fixed at the boundaries. The question is, however, whether the boundary terms 
in the fermion action need to be given "special" coefficients to avoid boundary errors 
(as we have seen for the gauge action). It turns out that this is not the case for the Wi 
action [J35[|. We will briefly discuss this issue later. For the moment, we will proceed by 
assuming that the Wi and W2 actions with Dirichlet boundary conditions defined in the 
previous paragraph only have the 0(a n ) on-shell errors of the bulk actions (with n = 2 
for Wi, n = 3 for W 2 ). 

One of the important features of the Schrodinger functional for practical applications 
is that with Dirichlet boundary conditions there are no zero modes at vanishing quark 
mass, at least classically. This is true even in the continuum limit with homogeneous 
Dirichlet boundary conditions, p = p = p' = p' = 0, where, at tree level, the lowest 
eigenvalue of the Dirac operator is ir/2T (one actually has to consider eigenvalues of 
75-^); see [19, 36] for details). 



3.3 Boundary Fields 

It is very useful (e.g. in the non-perturbative determination of the clover coefficient, 
see sect. |j) to define boundary fields || as functional derivatives with respect to the 
boundary values of the quark fields. 



19 



For the case of a Wi action these functional derivatives will always be with respect 
to the boundary values in eq. ( p. 12 ). To see explicitly what these boundary fields look 
like, consider the terms in a Wi action coupling to the (inner) boundary at xq = 0, say, 



« 3 E 



p(x)ao^p(x) - p(x)C/ (0,x)p-V(a ,x) - ^(a , x)P+C/ (0, x)tp(x) 



(3.13) 



There is a similar formula for the terms coupling to the boundary at xq = T. The 
correlation functions of the boundary fields we will consider later only involve the zero 
momentum Fourier components of the boundary fields. The spatial derivative term 
in ( |3.13 ) does not contribute then, and we will from now on drop such terms. 
Let us define the boundary fields by applying 

C(x) = T^-r , C(x) = 



C'(x) = 



<5p(x) 
5 

WW) 



C'(x) 



<5p(x) 
5 

WW 



(3.14) 



within the functional integral (as left-derivatives, say) and then setting all p and p to 
zero. Since not all components of the boundary values are independent (e.g. P p = 
from eq. fl3.12|) ), the above definition is ambiguous. To achieve uniqueness we impose 
the natural constraints 



P+C(x) 
P-C'(x) 



o, c(x)p- 

0, C'(x)P H 








(3.15) 



Explicitly, we then find 



C(x) 
C(x) 
C'(x) 
C'(x) 



E/b(0,x)p-V(ao,x) 
^(ao,x) P + C/ o (0,x)t 

c/ (r , -a ,x)tp+v;(r / -ao,x) 

4>(T'-a , X )p-U (T'-a ,x) . 



(3.16) 



For a W2 action one can in principle define two sets of boundary fields, either by 
taking derivatives with respect to the inner, or with respect to the outer boundary 
values. However, the latter are much simpler: The terms in the action coupling to the 
outer boundaries are given by (ignoring spatial derivative terms) 

b a 3 J] [p(x) U (0, x) P~ VK x) + ij(a , x) P+ U (0, x)* p(x) 

X 

+ /? , (x)C/o(T / -ao,x)tp+V;(T / -ao,x) + ^(T , -a ,x) P~ C/ (T'-a ,x) ^(x)] ,(3.17) 

where we have used the fact that the temporal links within the double boundary layers 
are all 1. We see that if we now define, with respect to the outer boundary values, 

5 ~z * \ 5 



C(x) = 
C'(x) = 



b 5p(x) 
S 

WW ' 



C(x) = 
C'(x) = 



b Sp(x) 
5 

WW 



(3.18) 
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the explicit expressions for the boundary fields are identical to those for Wi actions, 
eq. flOj)! 



3.4 Quark Propagator with Dirichlet Boundary Conditions 

The correlations functions we are interested in correspond to homogeneous Dirichlet 
boundary conditions, p = p = p' = p' = (remember that for the boundary fields we set 
the p's to zero after taking the functional derivatives). The propagator 

G(x,y) = mx)$(y)] F (3.19) 

in a given gauge background ([. . .}f denotes a fermionic contraction) therefore satisfies 

£ Q(x, z) ■ G(z, y) = 5 S xy , < x ,y < T' , (3.20) 

with the boundary conditions that G(x, y) vanishes when x or y are on either boundary. 



Eq. ( |3.20|) can be solved by standard iterative methods. 

It is easy to check that with homogeneous Dirichlet boundary conditions 75Q is 



hermitean (as long as M and Dk in eq. are hermitean, respectively, anti-hermitean) , 
so that 

G(x,yf = l5 G(y,x) l5 , (3.21) 
as for the usual case of periodic boundary conditions in all directions. 

We end this section with some brief remarks on the boundary errors of the fermionic 
part of the action. A detailed discussion will not be presented here, but we would like 
to point out the difference to the gluonic boundary errors. 

The main point is that Wilson-type quark actions are only on-shell improved, 
even classically and with periodic boundary conditions. This is due to the field 
redefinition necessary to avoid the doubler problem (cf. appendix C). One strategy for the 
determination of the coefficients of the boundary terms in the quark action is the study 
of correlations functions. Namely, to ask for which boundary coefficients the correlation 
functions of bulk fields, ip(x) and V>(x), as well as boundary fields, £(x) etc, are improved 
after appropriate field renormalizations have been performed. 

At tree- level it is sufficient to consider all possible two-point functions, by Wick's 
theorem. It is not difficult to see that with homogeneous boundary conditions the bulk 
field renormalization is given by the usual field redefinition factor, cf. appendix C. This is 
true for any value of the boundary coefficients in the action (since the bulk fields vanish 
at the boundary). The latter are constrained by the demand that the boundary fields 
renormalize multiplicatively. 

A detailed discussion of the two-point functions involving boundary fields can be 
found in |3E| for the Wi case. The conclusion is that the coefficients of the boundary 



terms in the action should be exactly as in the bulk (cf. subsect. 3.2); then the boundary 
fields are improved once they are renormalized by the inverse factor of the fields in the 
bulk. 
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For the W2 case the same is true for the zero momentum components of the boundary 
fields. When considering general boundary fields one has to appropriately adjust the 
coefficients of spatial derivative terms at the boundary (which we ignored, except in 
eq. ( 3.15| ), since they are irrelevant for the applications we have in mind). 

On the quantum level the general discussion of boundary terms at 0(a) is identical 
for Wi and W2 actions. We refer to H for the details. Of immediate interest to us is only 
the conclusion that the quantum boundary effects do not matter for the applications we 
are interested in. 



4 PCAC and Chiral Symmetry Restoration 

Here we sketch one of the simplest but most important application of the Schrodinger 
functional, namely the determination of the non-perturbative clover coefficient (s). 



Though much of the following is simply lifted from [14], we think it is at least 
pedagogically useful to go through this exercise, demonstrating, for example, how to 
maintain classical (plus tadpole, perhaps) improvement to higher order, while tuning 
quantum errors at 0(a). This includes some discussion of the anisotropic case. We will 
also provide a comparison of the Wi and W2 actions at tree-level. 

Consider lattice QCD with (at least) two flavors of mass-degenerate quarks. In the 
continuum limit we expect the isovector axial current to satisfy the PCAC relation 

8,Al = 2mP a , (4.1) 

with the associated pseudo-scalar density P a , as long as the quark mass m is sufficiently 
small. The idea behind using the PCAC relation in restoring chiral symmetry non- 
perturbatively, is the following. The PCAC relation is a local relation between operators 
(imagining, momentarily, that we continued to Minkowski space). It should therefore 
hold with the same mass m independent of what global boundary conditions are used,^] 
or what matrix elements of the PCAC relation are considered. As it turns out fln|, if 
one uses the wrong clover coefficient in the lattice action, m has a significant dependence 
on the kinematical parameters used to probe the theory. Requiring the dependence on 
kinematical parameters to vanish provides a powerful non-perturbative tuning method. 

To make these ideas more precise one has to consider the improvement and 
renormalization of the the axial current and density on the lattice (which can be 
anisotropic, at least for the moment). Naive, unimproved expressions for A^ and P a 
are given by 

P a (x) = V5(x) 75 ±tXx) . (4.2) 

Here r a is a Pauli matrix acting on the (suppressed) flavor indices of the quark fields J 1 ^ 

17 Note in particular, that the renormalization of the boundary terms in the gauge and quark actions has 
no effect [H], as already remarked in sect. 2. One can therefore use the classically improved Schrodinger 
functional when determining the clover coefficient non-perturbatively. 

18 We normalize the Pauli matrices such that Tr T a T b = 28 a h. 
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As shown in appendix C, on the classical level improvement of the operators A®(x) 
and P a {x) changes them only by the same overall mass-dependent factor, to the order 
we are working. The PCAC relation is therefore not affected. 

On the quantum level there will be additional errors starting at 0(a). One can 
show that the only modification of the PCAC relation required at 0(a) is, in continuum 
notation, the addition of a counter term 

(<L4g , SA%) = (o c A0 d P a , ac A d k P a ) (4.3) 

to the axial current A°(x)0 (see J5j for a thorough discussion of the renormalization of the 
PCAC relation). Here we are ignoring the multiplicative renormalization of A^ + SA^ 
and P a on the quantum level. For the non-perturbative determination of the clover 
coefficient the precise values of these renormalization factors are not required. We will 
absorb them in the unrenormalized current quark mass m (and continue calling it m). 

Consider now correlation functions of the PCAC relation with some product of fields 
O. Then 

V M ((A»+oA«)(x)0) = 2m(P a (x)0) + 0(a 2 ) , (4.4) 

if we have chosen the coefficients of the improvement terms equal to their non- 
perturbative values (for the W2 case the lattice derivative in this and other equations 
to follow is understood to be the improved one; cf. the last footnote). Note that O does 
not have to be a renormalized operator; the correction terms will cancel between the two 



sides of eq. (4.4). What is important, however, is that the support of O does not contain 
x, to avoid contact terms. Finally, we should remark that for the W2 case the a 2 errors 
in ( [4.4[) are pure quantum errors (for suitably defined m, cf. above). 

A good choice of O is given by a product of zero momentum quark boundary 
fields |, ||. Define, for X a = Ag and X a = P a , 

fx(xo) = ~a e Yl l(X a (x)((y)^r a az)), (4.5) 

in terms of the lower boundary fields £ and £. Due to spatial translation invariance the 
lhs depends only on x > as anticipated by our notation (in practice one should average 
over spatial x on the rhs to increase statistics). We now see that the PCAC relation 
implies for the current quark mass, 

m(x ) = r(x ) + a CAos(x ) = m + 0(a 2 ) , (4.6) 

where 

K*„) S |^ff .(„) S (4.7) 

Note that due to spatial translation invariance the dependence on ca has dropped out 
of eq. (E 
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In the W2 case an improved discretization of the derivative in (4.3) has to be used on the lattice, 



which, borrowing notation from covariant lattice derivatives, we can write as V M (1 — ^alA 
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So far we have considered the general case of anisotropic lattices. The isotropic 
case is significantly simpler in practice, since then only uo = ojq has to be tuned in the 
quark action. In the anisotropic case the color- electric background field chosen here 
should essentially couple only to the temporal clover term, ojq oofc FqHi an d allow its non- 
perturbative tuning, even if the spatial clover coefficient u is not exactly at its correct 
value. In a second step one could tune u by choosing boundary conditions corresponding 
to a color-magnetic background field (if necessary, one could then perform a second 
iteration of tuning ljq). 

With spatially periodic boundary conditions a background magnetic field is 
necessarily quantized with a minimum non-zero value that might well be too large 
(cf. e.g. |37]] and references therein) for our purposes. To avoid this problem one 
might want to choose suitable fixed boundary conditions in one (or more) of the spatial 
directions, which allows one to generate an arbitrarily small (classical) background color- 
magnetic field. In other words, one basically exchanges the time direction with one of 
the spatial directions in the Schrodinger functional. 

Note that anisotropic lattices have a distinct advantage in the context of Schrodinger 
functional simulations: One can use a small temporal extent T in physical units to give 
a large lowest eigenvalue (and thereby, presumably, alleviate problems with exceptional 
configurations), and still have a large number, T/clq, of time slices required from a 
practical point of view. 

These ideas should be subjected to some numerical studies (in the free and interacting 
cases). Since we have not done so yet, we will in the remainder of this section restrict 
ourselves to isotropic lattices. 

Let us continue describing the improvement conditions for determining the clover 



coefficient u. To eliminate ca = cao from (4.6) one can consider the analogously defined 



m'(xo), r'(xo) and s'(xq) in terms of the upper boundary fields £' and Namely, set 
A(T'-xo) = +a 6 £ I<^(zK'(y)75§r a C'(z)> 
f' P (T'-x ) = -a 6 E l(P a (x)C'(y) 1 ^rX(z)) (4.8) 

a,y.z 

and define r'(xo) and s'(xq) exactly as in ( [4.7D with /^(xo), fp(xo) replaced by f' A (xo), 
f'p(xo), respectively. The equations for m(xo) and m'(xo) now imply that 

CA{x ,y ) = —■ — — 4.9) 

a s(x ) - s'(y ) 

is an estimator (in the sense of probability theory) of ca, up to 0(a) errors. 

Before proceeding we have to specify the boundary conditions. For the Wi case 



a good choice of "boundary angles" (cf. ( 2.21 )) for the determination of the clover 
coefficient was found to be fb 



' ' ' (-hhh)* (4-10) 
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with vanishing quark "twist" , Ok = (cf . appendix A) . With these values the background 
field strength Fq^ is still small in lattice units, for typical L and T used in practice, but 
large enough to provide a significant lever arm for the tuning of uj. It is crucial that the 
boundary conditions are different at the upper and lower boundary; otherwise one could 
not use eq. ( |4.9| ) to eliminate c A from the problem of the uj determination. 



With the above boundary values it turns out that to keep the denominator in (4.9) 
as far away from zero as possible, it is advisable to use small xq and yo. They should 
not be too small, of course, to avoid boundary effects. After some study, the choice 
x o = Ho = \T was suggested [14] for the Wi case on an T • L 3 = 16 • 8 3 lattice. 



Let us now introduce the following observables as estimators of the (unrenormalized) 
quark mass, 

M(x ,y ) = r(x ) + ac A (yo,yo)s(x ) 

M'(x ,y ) = r'(x ) + ac A (yo,yo)s'(x ) . (4.11) 

One can then employ the demand that AM(xo,yo) = M(xo,yo) — M'(xo,yo) vanish, 
for suitably chosen xq and yo, as an improvement condition to determine uj |jl4| . Since 
AM(xo>^o) = 0, one should try to separate xq and yo a s much as possible. Using 
AM (IT, \T) is a good choice for the Wi case. 

We should mention a few more details: 

• The tuning of uj should ideally be performed at zero quark mass. Within the a 2 
ambiguity of the quark mass one could use any M(xo,yo) or M'(xo,yo) to define 
the quark mass. A good choice for the Wi case is M = M(hT, jT). 

• In practice it turns out that AM has only a very weak dependence on M. This is 
true in the free case for both the Wi and W2 actions up to quark masses as large 



as aM = 0.1 or 0.2. In [14] this was also found to hold for the interacting case 
(although only much smaller quark masses were investigated). One can therefore 
tune uj at small positive values of the quark mass, which alleviates problems with 
accidental quark zero modes on coarse lattices. It remains to be investigated exactly 
how large a mass one can use in the interacting case (cf. the end of this section for 
some relevant results). 

• AM does not vanish in finite volume, even at tree-level and zero pole mass. As 
improvement condition one should therefore not demand that AM vanish, but 
rather that it equals its tree-level value AM' ' . This is a relatively small correction, 
but should be taken into account, nevertheless, because it guarantees that tuning 
at positive quark masses does not lead to additional tree-level errors in uj (note 
that in such a case AM' ' has to be evaluated at a mass that is correct at least at 
tree- level) . 

• Due to the higher derivatives present in a W2 quark action one expects the 
boundary effects to extend one more lattice spacing into the bulk. For the W2 case 
one should therefore presumably move xq and yo in the improvement condition 
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AM(xo,yo) = A.M(°\xQ,yo) slightly further away from the boundary than in the 
Wi case. The necessity to do so can already be seen in the free case; the exact 
values to use for x$ and yo should be fixed after some numerical studies have been 
performed. 

To determine to one uses the estimator ( |4.9| ) to eliminate ca from the equations 
for m(xo) and m'(xo). This is adequate for the uj calculation, but for a more precise 
determination of ca itself it is better [[TJJ] to compare correlation functions with quark 
twist 6k 7^ with ones where 9^ = 0, in both cases with vanishing background field, 
4>k = <P'k = 0- More precisely, one should equate the corresponding mass difference to 
its tree-level value. It was found @] that this mass difference is again only weakly 
dependent on the quark mass (at least once the improved action with the previously 
determined clover coefficient is used). 

To give the reader some feeling for the various quantities involved in the improvement 
conditions for u> and ca we plot the tree-level current quark masses with the relevant 
boundary conditions in figures [| and ^, for both the Wi and W2 actions. Recall that at 
tree-level ca = 0, so that m(xo) = r(xo) and m' (xo) = r' (xo) . 

We should comment on one aspect of figure p|. Note that for the W2 action with 
u = the current quark mass estimates t(xq) and r'(xo) deviate less from the pole mass 
than for the corresponding Wi case. The reason for this is simply the smaller value 
of the Wilson parameter (r = 2/3 versus r = 1), making the W2 action somewhat less 
sensitive to errors involving the background field. This suggests to use a slightly larger 
background field when tuning a W2 action. 

For completeness we end this section by describing how to calculate the correlation 
functions /a(^o) an d fp(xo) in practice (analogous remarks hold for f A (xo) and fp(xo)). 

Writing the operator X in the form 'tp(x) ^T a Fx ip{x), we have 

fx(x ) = -a 6 E i(^(x)Ir a r x V(x)C(y)75^C(z)) 

= +|a 6 ^ (Tr 75 [C(z)^(x)] F r x ^(x)C(y)]F>G J (4.12) 

where (. . .)q denotes the average over gauge configurations (the trace, to be sure, is over 
spin and color). From (|3.16j ) we read off 

[C(z)V(x)] F = Uo(P,z)P-G(z 

) x )\z =a , 

mx)C(y)} F = G(x,y)\ yo=ao P+U (0,y)K (4.13) 
let us define 

H( x ) = a 3 ^G(x,y)| yo=ao P+[/ (0,y)t, (4.14) 
y 

which is the solution of 

Y d Q{x,z)-H{z) = -5, , ao P+C/ (0,x)t, 0<x <T', (4.15) 
a 



To simplify (4.12 
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Figure 5: The dimensionless ratios r(xo)/m p (left) and r'(xo)/m p (right) as a function of 
xq/cl, for the free SW (or W±, +) and D234 (or W2, •) actions on an isotropic lattice. m p 
is the pole mass, the zero momentum energy of a free lattice quark in in the absence of a 
background field. The points in the upper/lower half of the left/right figure correspond 
to vanishing clover coefficient uj = 0; the others to classical improvement, u = 1. The 
kinematical parameters used are, for both actions, am c = 0.02, T = 16a, L = 8a, 6^ = 0, 
with the angles 4>, cf)' specifying the background field given by eq. ( [4.10 ). 
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Figure 6: As in figure || but with vanishing background field, (p = <fr' = 0. Now the 
points in the upper half correspond to 0k =ir/3, the rest to 8k = (in the absence of a 
background field r(xo) = r'(xo), and the value of uj does not matter at tree-level). 
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and the boundary condition that H (x) vanishes on the boundary layers. Physically H(x) 
is the propagator for a quark from anywhere on the lower (outer) boundary to the point 
x in the bulk. Note that HP~ = 0, so only half of the components have to be calculated. 
Using ( |3 . 2 1| ) one finally obtains 

f x (x ) = i(Trtft( x ) 75 r x #Or)} G . (4.16) 

Note that j§Tx = 1 for X = P, which means that fp(xo) is positive configuration by 
configuration. 



We used these formulas for the free case, solving (4.15) with standard inversion 
algorithms, to obtain the numbers plotted in figures I andg 

As a check of some of our earlier remarks (and to see how small a volume can be 
used in practice), we have repeated the determination of the non-perturbative clover 
coefficient of the SW action on quenched (5 = 6.0 Wilson glue on a 12 • 6 3 lattice at 
aM = 0.04. We obtain uj = 1.74(4), in agreement with the original calculation on a 16 -8 3 



lattice at aM = 0.01 [14 1 



5 Conclusions and Outlook 

We are currently witnessing dramatic progress in lattice gauge theory, related to our 
increased understanding of the design of improved actions. More precisely, there are two 
important developments. One is the ability to perform simulations on coarse lattices of 
improved glue, the other the non-perturbative elimination of O(a) errors in Wilson-type 
quark actions (and operators), as suggested by the ALPHA collaboration. Unfortunately, 
with quenched Wilson glue the non-perturbative improvement of the SW action is not 
possible on lattices coarser than 0.1 fm; the quenched approximation seems to break 
down for small quark masses on such lattices. One thereby looses the great advantage 
of coarse lattices, the possibility of performing a series of relatively cheap simulations 
on such lattices that can then be accurately extrapolated to the continuum (the latter 
assumes, of course, that O(a) errors are negligible). 

The aim of this work has been to set the stage for the use of non-perturbatively 
improved quark actions with improved glue, by adapting the Schrodinger functional to 
such actions. Much remains to be done. We have sketched only one application, the 
tuning of the clover coefficient(s), and even this area requires further study, certainly on 
anisotropic lattices. One-loop calculations, although more difficult for improved gauge 
actions, should eventually be performed for various improvement coefficients. 

The motivation for this work has been the hope that it is possible to non- 
perturbatively improve quark actions on significantly coarser lattices when using 
improved rather than Wilson glue. This indeed seems to be the CclSG, clS will be 
described in |25| for the case of the SW action with a 2 improved glue (in the quenched 
approximation). Besides the SW action one should also consider the D234 quark action, 
which has much smaller violations of rotational symmetry. It will be interesting to 
see to what extent the occurrence of "exceptional configurations" is affected by the 
choice of quark action. As mentioned in the introduction, our own experience with 
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spectroscopic simulations indicates that the D234 action has less near-zero modes (on 
the same configurations, at the same quark mass). 

Another option is the use of different kinds of improved glue. For the Wilson gauge 
action the breakdown of the quenched approximation occurs on lattices that are still 
surprisingly fine (at least it seems so to us). A detailed understanding of this breakdown 
does not exist yet, but it might, at least partially, be related to the (non-perturbative) 
phase structure of the Wilson action. If such strong-coupling details of the action 
are relevant, certain types of improved glue might allow non-perturbative tuning on 
significantly coarser lattices than others. 

The problem of "fake" near-zero modes does not exist in the continuum limit; their 
occurrence is a lattice artifact related to the specific Wilson-type discretization used 
for the quark action (the problem does not occur for staggered fermions). It would 
therefore not really be surprising if improved actions have less problems with exceptional 
configurations; such actions are, after all, supposed to be more continuum-like. Note, for 
example, that adding higher derivatives to an action makes the typical paths contributing 
to its path integral smoother (cf. e.g. |38|]).0 Of course, this statement is rather formal, 
and numerical studies of these questions are called for. 

A different strategy is to try to modify the quenched approximation to eliminate 
the problem of exceptional configurations. The challenge is to find a method that is 
systematic, practical, and satisfactory on a conceptual level. For a specific suggestion 
and some encouraging numerical studies, see |Q. Such a procedure will most likely 
allow non-perturbative tuning on much coarser lattices than possible within the standard 
quenched approximation. 

We have discussed the 0(a) improvement of quark actions on anisotropic lattices, and 
referred to the decisive advantages such lattices have for the study of heavy quarks^ and 
glueballs, for example. We suggested how to use the Schrodinger functional to achieve 
0(a) improvement also on anisotropic lattices. 

Finally, due to the great cost of full QCD simulations we regard it as very important 
to study the improvement of quark actions coupled to dynamical improved glue, on 
isotropic and eventually also anisotropic lattices. 

We hope that the program outlined here will lead to the construction of actions that 
allow accurate simulations of all aspects of QCD on coarse lattices (this will also require 
the determination of various operator renormalization constants along the lines pioneered 
by the ALPHA collaboration). The advantages of coarse lattices make it imperative, we 
think, to pursue these studies. 

20 Improvement of a Wilson-type action at O(a) is an exception in this respect; adding the clover term 
leads to more, not less, problems with exceptional configurations [Q. It is only higher derivatives in the 
kinetic term that lead to smoother typical paths. 

21 Note that there are signs E9| that NRQCD breaks down for charmonium. The D234 action on 
anisotropic lattices would appear to be an ideal tool || [lJJ for the study of charmonium and heavy-light 
systems. 
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Appendix A Notation and Conventions 

A.l In the Continuum 

We write the action of euclidean SU(iV) gauge theory in four dimensions as 

S g [A] = ^Jd 4 x Tr F^(x)F^{x) , (A.l) 

where F^ u (x) is the su(iV)-valued hermitean field strength. In our conventions the 
covariant derivative is expressed in terms of the hermitean gauge field A fJi (x) as 

D„ = d^-iA^ (A.2) 

so that 

Ffj, v = i[D^D v ] = d il A v -d v A il -i\A tl ,A v \ . (A.3) 

We use traceless hermitean su(N) generators T a , a = 1,...,N 2 — 1, normalized by 
Tr T a T b = \5 ab . We write F^{x) = F* v (x)T a and A^x) = A*(x)T a . 
The action of a Dirac fermion coupled to a gauge field is 

S q $, i>,A] = J d A x 4>{x) {p + m) tf>(x) , (A.4) 

where Ip = Ylu 7m-^V * n terms °f the euclidean gamma matrices defined by 

7^ = 7j 

{ 7m , 7 ,} = 25^ . (A.5) 

The spinor fields ip(x) and 4>(x) also carry a suppressed color index in the vector 
representation of SU(AQ, and perhaps another suppressed flavor index. 

The hermitean matrices are defined by a^ u = — §[7^,7^] or 7^7^ = 5^ u + ia^. 
Note that 

Ip 2 = J; Dj + i^a^. (A.6) 

(J, /IV 

The term a-F = cr^F^ is known as the clover term (cf. below). 



A.2 On the Lattice 

We will consider a four-dimensional euclidean hypercubic lattice of extent and lattice 
spacing m direction (j, = 0,1,2,3.0 O n a spatially symmetric lattice L = for 
k = 1,2,3. For purely classical considerations the are arbitrary, but for any situation 
where quantum effects play a role we will assume that the lattice is spatially isotropic, 
a = a/-, k = 1, 2, 3. 

22 Note that depending on how one codes the Schrodinger functional the number of time slices used in 
one's program, Lo, will differ from the "true physical extent" T of sect. ti. 
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Points are labelled by x,y,z .. ., as in the continuum. Spatial vectors are denoted by 
boldface letters, as in x = (xo, x). We use k, I for spatial indices. The notation x ± is 
a shorthand for x ± a^fi, where /t is a unit vector in the positive //-direction. 

When working on anisotropic lattices one should in principle be careful in specifying 
the lattice spacing errors of various quantities. To avoid cumbersome notation, we will 
use 0(a n ) to denote any combination of errors of the form 0(a"), up to logarithmic 
corrections. Only occasionally do we use the notation 0(a") for an error depending only 
on a M , for some specific fi. 

The lattice gauge field, or link field, will be denoted by U^x) or U^{xq,sC). 
Mathematically it is a parallel transporter from x + fi to x, along a straight line (the 
link). In terms of an underlying continuum gauge field Aa(x) it can be written as 



UJx) = Pexp 



dx' A^x') 



(A.7) 



where P denotes path-ordering, stipulating that fields later on the path are to be placed 
to the right of earlier fields. We will employ the notation U-^(x) = U^(x — /x)t for the 
parallel transporter from x — fi to x. 

Our notation for lattice covariant derivatives was presented in sect. 3. We should 
add one detail here ||. If one wants to impose "twisted" boundary conditions on the 
quark fields, 

ip{x + L k k) = e i9k ip(x), i)(x + L k k) = e~ i9k ip(x), (A.8) 

one can, equivalently, multiply each link field Uk(x) in the derivative operators of sect. || 
by a phase exp^a^^/L^), while keeping the quark fields strictly periodic. This twist 
corresponds to adding an abelian gauge field to the action, so that, for example, the 
covariant continuum derivative is given by = <9 M — iA^ + iO^/L^ (6>o = 0). 

When working up to 0(a 3 ) errors in the quark action, classically, we can use the 
clover representation of the lattice field strength (see e.g. Q) in the a-F term. This 
discretization differs at order a 2 from the continuum field strength. To eliminate the 
a 2 errors in the field strength one can use an improved discretization, for which we also 
refer to B. 



Appendix B Classical Boundary Errors of the Pure 

Gauge Schrodinger Functional 

We want to study the difference between the lattice and the continuum gauge actions 
with corresponding boundary conditions. As remarked in sect. this difference has two 
sources; one being the difference between an integral and a sum, the other the difference 
between the integrand and the summand. 

When considering the continuum limit, classically or in perturbation theory, one 
assumes the lattice field to be given by U^x) = Pexp[— i J^ +fl dx'A^x')], in terms of a 
smooth continuum field A^(x) with the appropriate boundary conditions. Recall that the 
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difference between the integral and a sum over a smooth function with periodic boundary 
conditions is exponentially suppressed in the lattice spacing; the exponential decay rate 
being equal to the minimal distance of the singularities of the integrand from the real 
axis (this follows from the Poisson summation formula). 

In the bulk, the lattice action densities differ from the continuum one at the order 
a n of improvement. Any larger errors will be due the difference between the sum and 
an integral at the boundary, or, possibly, due to a difference between the continuum and 
lattice integrands at the boundary (that is larger by more than a power of a than the 
difference in the bulk, since the effects of a boundary of codimension 1 are suppressed 
by such a power). We will see that for improved actions both errors contribute — and 
in fact partially cancel. 

It will be useful to first describe the simpler case of the Wilson action. In this case the 
difference between the lattice action density and the continuum density cx Tr F^ u (x) is 
0(a 2 ) everywhere, including the boundary, if we always associate the Tr Fj% u (x) arising 
from the expansion of a plaquette to the middle of the plaquette. For the temporal 
integraiP| over the F kl terms arising from spatial plaquettes within time slices, we then 
have a sum of the general form (cf. figure |2[) 



hf(0) + f(a) + ... + f(T-a) + hf(T) 



T 



dtf(t) + 0(a 2 



(B.l) 



(to avoid cumbersome notation we use a for the lattice spacing, even though it should 
more properly be do for a temporal integral). This is the trapezoidal rule. 

Similarly, for the F^ k terms arising from temporal plaquettes living between time 
slices, the midpoint rule tells us 

a [/(a/2) + /(3o/2) + . . . + /(T-3o/2) + /(T-a/2) 



f T dtf(t) + 0(a 2 ). (B.2) 
Jo 



All in all we have shown that for the Wilson case there are classically only 0(a 2 ) boundary 
errors, if we choose the temporal weight factors as in figure [2|. This was first shown in (T^j 
by a different method. 

Before proceeding with the improved case let us point out that the sum rules flB.l|) 
and (|B.2j) can be derived by "stringing together" their elementary building blocks 



hf(0) + \f{a) 



dtf(t) + 0(a 3 ) 



(B.3) 



and 



a /(a/2) = f a dtf(t) + 0(a 3 ), 
Jo 



(B.4) 



respectively. These formulas are of course trivial; they and their slightly more 
sophisticated cousins used below can all be derived in the same manner: Since the 



Due to the spatially periodic boundary conditions the space integral always comes along for the ride 
and can be ignored. 
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integrand can (by assumption) be taylor expanded, one only has to prove them for 
polynomials up to the appropriate order. 

It will probably not have escaped the reader's attention that there is a complete 
analogy between the elementary sum rules and their extended forms on the one hand, 
and the transfer matrix for one time step and T/a steps on the other. 

Let us now apply the above lesson to the improved case. First of all note that the 
continuum boundary layers should be placed in the middle of the boundary double layers 
of figure We will relabel the time axis in this figure so that the continuum boundaries 
are positioned at t = and t = T. The lattice time slices now run from t = —ao/2 to 
t = T + in steps of size ao (which will be denoted by a again, from now on). 

Consider first the errors associated with the spatial loops. Up to 0(a 4 ) errors we 
can associate an F% 1 term to the middle of each spatial plaquette (by construction of 
the improved action). As in the Wilson case, we can now proceed to analyze the sum 
over time slices for fixed spatial x. The question is whether there exists an appropriate 
midpoint rule approximating an integral up to 0(a ) errors. The answer is positive, 

1 23 23 1 

24/(-a/2) + -/(a/2) + /(3a/2) + . . . + f(T-3a/2) + -f(T-a/2) + -f(T+a/2) 



[ T dtf{t) + 0(a 4 
Jo 



(B.5) 



which follows from the elementary rule 



l/(-a/2) + §/(a/2) + ^/(3a/2) 



f a dtf(t) + 0(a 5 ). (B.6) 
Jo 



Note that it is a priori quite surprising that it is possible to integrate the polynomials 
l,t, t 2 ,£ 3 exactly with only three sampling points. In any case, we have shown that one 
should choose b = 1/24 in figure 

We now come to the errors associated with temporal loops. These are a bit more 
subtle. The reason is the following: In the bulk the plaquettes and rectangles combine 
to give, up to 0(a 4 ) errors, an F^ k term situated at the middle of each plaquette; just as 
for the spatial loops. However, this not true for the boundaries Q which introduces an 
a 2 error. Note that this error is associated with the integrand at the boundary. There 
is another a 2 boundary error, due to the difference between the trapezoidal rule and the 
continuum integral. It turns out that the corresponding correction term, exhibited in 
the improved trapezoidal rule, 



5/(0)+/(a) + ...+/(T-o)+i/(r) 



12 



(/'(T)-/'(0)) 



l T dtf(t) + 0(a 4 ) (B.7) 
J o 



is exactly cancelled by the error from the integrand at the boundary! We leave the explicit 
proof of this statement as an exercise to the reader (who might want to consult || or |41], 
for example, for the small a expansion of the plaquette and rectangle terms). 

24 It were true if each of the tall-temporal rectangles that overlaps with the boundaries in figure ^ were 
split into two copies (each copy given a weight of |), and one copy shifted by one temporal lattice spacing 
so that it sticks out beyond the outer boundary layer of the lattice. 
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For completeness we remark that the improved trapezoidal rule follows from the 
elementary rule 



1 13 13 1 



[° dtf{t) + 0(a 5 ) (B.8) 
Jo 



(or from the Euler-Maclaurin formula). This completes our proof that with b = 1/24 in 
figure |2] the Schrodinger functional for improved glue has no classical errors larger than 
0(a 4 ). 



Appendix C Quark Actions on Anisotropic Lattices 

Consider an anisotropic lattice that is spatially isotropic. Our main aim in this appendix 
is to show that an on-shell 0(a) improved Wilson- type quark action on such a lattice 
can always be written as a discretization of an effective continuum action of the formP] 



S = d 4 x ip(x) m + IpQ + cp - \ra Q ( ^ Df t + u> a ok Fok + oj a k iF M j 
J L m k k<i ' . 

(CI) 

Here is the covariant (continuum) derivative, and c, ujq and u are the parameters that 
have to be tuned at the quantum level, c is a bare velocity of light that has to be adjusted 
so that the renormalized velocity of light is 1, e.g. by studying the pion dispersion relation 
at small masses and momenta, ooq and u> are the temporal, respectively, spatial clover 
coefficients. Classically all three coefficients are 1, independent of the value of the Wilson 
parameter r. On an isotropic lattice c = 1, and only ujq = to has to be tuned. 

To avoid doublers, the Wilson parameter can in principle be chosen to have any 
value r > (for the W\ case there is an upper bound of 0(1) from requiring reflection 
positivity; the exact value depends on the specific discretization used). Conceptual and 
practical considerations usually lead to a specific choice of r for any given discretization; 
cf. (|3TT0|) . The non-perturbative values of ujq, uj and c depend on r, of course. 



Note that to 0(a) we can without loss of generality always assume that ujq and to are 
mass-independent. On the other hand, c can have a mass dependence on the quantum 
level at 0(a). 

Before proving the above claims, it will be useful to provide a bit of background on 
classically improved, doubler-free quark actions on anisotropic lattices. The simplest way 
to construct such actions || is to start with a continuum action ip c Q c tp c = ip c {]p + m c )Tp c 
and perform the field redefinition vp c = Q c i/j, V> c = i/j Cl c with 

n c = n c , n c n c = 1 - ^ (p-m c ) . (02) 

The subscripts 'c' stand for 'continuum'. A field redefinition is just a change of variable 
in the path integral, so on-shell quantities are not affected .p 5 ] The transformed quark 

25 A similar statement, in a different but refated context, can be found in 

26 The Jacobian of a field transformation matters only at the quantum level, where, in the case at hand, 
its leading effect is to renormalize the gauge coupling. 
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matrix reads 



(l c Q c n c = m c (l + \ra m c ) + Jf> - - ra D 2 ^ + \a-F) , (C.3) 



where we used eq. (p 

Note that the above still defines a continuum action; ao should at this point just 
be considered a formal dimensionful parameter. It takes on the meaning of a temporal 
lattice spacing only in the next step: We obtain the Wi and W2 actions of sect. || simply 
by replacing the continuum operators Ip, clqD 2 and a^F^ by suitable lattice versions, 
differing at 0(0") from the former, with n = 2 for Wi and n = 3 for W2. For a more 
detailed exposition see ||. 

By construction, classical lattice actions obtained in this manner have on-shell errors 
only at 0(a n )J^| The proves our claim relating to (C.l) in the classical case. 



Note that 0(a n ) errors can also be achieved for off-shell quantities simply by undoing 
the field redefinition (on the lattice now), to obtain continuum- like lattice fields. By abuse 
of notation we will also denote the lattice fields for the action obtained by discretization 
of ( |C.3| ) as i/j, (like the fields defined by the inital change of variable in the continuum). 
For quark bilinears, where things are particularly simple, we can then rewrite bilinears 
of the lattice fields in terms of continuum fields as 

tp(x)Fip(x) = $ c {x)q.- 1 tq.- 1 $ c (x) + o(a n ) 

= {l + ra G m c )- l ^ c {x)Tij c {x) + 0{a n ) . (C.4) 

In the above we used the field equation for ip c to eliminate lp in the field transformation 
operators. We see that for bilinears classical improvement simply amounts to the 
multiplication by a mass-dependent factor (that, furthermore, is the same for all 
bilinears). 

Let us now return to our claim for the quantum case at 0(a). The possible terms 
that can appear at this order in the effective action are of the form 

Dl J2 D l W>>Pol J2 a okF ok , J2 a kiFki- (C5) 

k k k<l 

Remember also that at 0(a°) the relative coefficient of the temporal and spatial 
derivatives can renormalize. 

To prove our claim it suffices to show that starting with the continuum action 
/ Tp{m c + Pq + c'~tf)]ip we can find a field redefinition that generates the D 2 , D 2 and 
[rJ,^o] terms in ( |C.1| ) with arbitrary coefficients; then we can always adjust the values 
of these coefficients to equal those in ( |C. 1| ) (e.g. zero for the [y>,flo] term). In other 
words, Dq, D\ and [r^,^o] are redundant operators. 



27 Note, in particular, that if one writes mo = m c (l + ^ram c ) in eq. (3.2), then m c agrees with the 
pole mass up to a n errors. 
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We assume that we have used our freedom of rescaling the fields to ensure that the 
coefficient of flo is 1 at 0(a°). The most general field redefinition at 0(a) is then 



£7 = 1 — 



n = 1 - 




r m ao 
4 



r m ao 
4 



(C.6) 



This leads to 



tt (m c + po + c'p) tt 

= m c [l + ^-m c (r m + f m )} + p [1 + ^-m c {r m +f m -r - f )] 



+ ^[c' + -m c (c'(r m + f m ) -r-f)] - -j(r + f )^ -- 'y 

- ^ [c'a (r + f ) + o(r + f)] # } 

- J [c'a (r - f ) - a(r - f )] [0, # ] ■ 




(C.7) 



Remembering J/) 2 = T,k D l + Y.kKi^kiFki and {^,$>o} = EfcCofc-^ofc, we see that our 
claim follows if we parameterize the independent free parameters of the change of variable 
as ro + fo (to adjust the Dq coefficient), r + f (to adjust J2k^l), (*° adjust [p,$>o]) 
and r m + r m (to eliminate any possible mass dependence of the po coefficient at 0(a)). 
No further simplifications are possible by suitable choice of tq — fo or r m — f m ; the latter 
does not even appear in the above. 

The mass dependence of c at 0(a) is a complication in practice. To minimize 
systematic uncertainties it is preferable to always determine the clover coefficients ujq 
and oj at zero quark mass. To do so one needs to know c at zero quark mass. It remains 
to be seen if the most naive method of determining this value — demanding the pion 
dispersion relation to be relativistic at small momenta, and then extrapolating to zero 
pion mass — is sufficiently accurate. One expects that c has only a weak dependence 
on ujq and uj; otherwise, one might have to repeat the above procedure iteratively until 
self-consistent values of c, ujq and uj are obtained. 
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